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摘 要 在 评估 国际 常用 Kp 指数 现 报 模式 Takahashi 及 其 应 用 于 中 国 地 磁 台 站 效果 的 基础 上 , 提出 了 一 种 改进 
的 地 磁 Kp 指数 现 报 模式 , 其 可 以 有 效 识 别 地 磁 规 则 日 变化 的 逐日 变化 特性 , 反映 地 磁 扰 动 的 季节 效应 和 地 方 时 效 
应 , 从 而 提升 了 Kp 指数 现 报 的 准确 性 . 采用 Takahashi 模式 开发 时 所 使 用 的 台 站 数据 进行 对 比 , 新 模式 将 单 站 地 
磁 Kp 指数 现 报效 率 由 0.77 提升 至 0.84, 多 站 联合 Kp 指数 现 报效 率 由 0.88 提升 至 0.92; 采用 2000—2006 年 北 
京 十 三 陵 台 站 (BMT, 磁 纬 29"N) 数据 进行 评估 ，Kp 指数 现 报效 率 由 0.70 提升 至 0.80. 对 Kp 指数 现 报 结果 的 
误差 分 析 发 现 , 现 报 误差 存在 明显 的 地 方 时 差异 和 一 定 的 季节 差异 , 误差 随 扰动 强度 变化 并 在 中 强 磁 扰 时 最 大 . A 
用 SuperMAG 的 指数 分 析 表 明 ，Kp 台 站 的 经 度 不 均匀 分 布 会 对 现 报效 果 造 成 一 定 影响 . 
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Abstract The planetary three-hour-range Kp index is widely used in space weather services. Due 
to the two weeks’ time delay of the release of Kp index, the real-time estimation of the Kp index 
becomes essential for space weather forecasting organizations. In this paper, an improved algorithm 
for real-time Kp estimation on the basis of Takahashi's procedure was developed. The new algorithm 
improved the accuracy of Kp estimation by distinguishing the day-to-day variations of geomagnetic 
field's regular variation and taking account of both the diurnal and seasonal variations of magnetic 
disturbance. The statistical result shows that the Prediction Efficiency (PE) increased from 0.77 to 
0.84 using Fredericksburg's data, from 0.88 to 0.92 using the nine stations network's data. Using the 
data from Beijing Ming Tombs observatory (BMT) from 2000 to 2006, the PE increased from 0.70 
to 0.80 compared with the Takahashi procedure. The error of the estimations differs from universal 


time, season and the scale of magnetic disturbance. The RMS error reached maximum for Kp = 7. 
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Analysis of Kp using the SuperMAG indices indicated that the asymmetric distribution of Kp stations 


affects the Kp nowcast performance. 


Key words 


0 引言 


国际 3h 磁 情 指数 Kpl- 是 磁 层 物理 研究 的 重 
要 参量 4-12, 同时 , 由 于 地 磁 扰 动 和 其 他 磁 层 物理 
现象 的 紧密 关系 , 很 多 空间 天 气 模式 需要 以 Kp 指数 
作为 输入 参数 ， 学术 研究 和 空间 天 气 预 报 对 Kp fü 
数 发 布 的 时 效 性 要 求 是 完全 不 同 的 . 目前 , Kp 指数 
的 发 布 (德国 波茨坦 地 球 科学 研究 中 心 ) 延迟 2 Jj, 
这 一 延迟 对 空间 天 和 气 预报 的 实时 性 要 求 影响 很 大 . 

当前 解决 Kp 指数 在 空间 天 气 预报 中 时 效 性 需 
求 的 方法 大 致 归 为 两 类 . 一 是 建立 Kp 指数 短期 预报 
模式 , 以 ACE 卫星 工 1 点 的 太阳 风 参 数 为 输入 , 提前 
约 半 小 时 预报 地 磁 Kp 指数 . 此 类 模式 主要 是 基于 神 
经 网 络 方法 开发 的 , 例如 Wing 等 的 APL 模式 U91, 
Boberg 模式 [和 Bala 模式 中. 二 是 根据 地 磁 台 站 
实时 监测 数据 发 展 的 地 磁 Kp 指数 佑 计 方 法 , 对 Kp 
指数 进行 现 报 . 此 类 模式 一 般 是 按照 Kp 指数 的 推导 
过 程 , 在 某 些 环节 上 采用 一 些 近似 假设 , 进行 Kp 指 
数 的 估算 , 具有 代表 性 的 是 Takahashi 等 16 采用 9 
个 地 磁 台 站 的 历史 数据 构造 的 现 报 模式 , 单 站 Kp 指 
数 现 报 值 与 实测 值 的 相关 系数 为 0.85~0.9, 而 9 台 
站 联合 Kp 指数 现 报 值 与 实测 值 的 相关 系数 为 0.94. 
该 模式 在 美国 国家 海洋 和 大 气管 理 局 空间 天 气 预 警 
中 心 (NOAA/SWPO) 运行 , 并 作为 地 磁 暴 警报 的 判 
断 依据 . 

Kp 指数 现 报 的 误差 来 源 主 要 有 两 方面 : 一 是 现 
报 模式 中 近似 的 假设 , 主要 涉及 地 磁 规 则 日 变化 提取 
以 及 由 单 站 或 多 站 天 指数 估算 Kp 指数 时 的 拟 合 关 
系 ; 二 是 用 于 Kp 指数 现 报 的 台 站 数 及 地 理 分 布 . 

降低 现 报 模式 误差 , 提升 地 磁 Kp 指数 现 报 的 准 
确 性 是 空间 环境 预报 的 重要 需求 ， 由 于 中 国 台 站 磁 
纬 较 低 , 地 磁 扰 动 的 响应 相对 高 纬 台 站 较 弱 , 当 完 全 
使 用 Takahashi 模式 时 , 正午 附近 现 报 值 与 Kp 的 相 
关系 数 较 低 . 此 外 ，Takahashi 模式 中 的 一 些 近 似 假 
设 使 得 该 模式 无 法 去 除 地 磁 规 则 日 变化 的 逐日 变化 ， 
对 地 磁 扰动 的 季节 变化 也 未 加 区 分 , 这 都 给 该 模式 的 


Kp index, Są index, FMI method, Nowcast, Magnetic disturbance 


地 磁 Kp 指数 现 报 准确 性 带 来 影响 . 

本 文通 过 对 以 Takahashi 模式 为 代表 的 国际 常 
用 Kp 指数 现 报 模式 计算 流程 的 分 析 , 对 该 模式 应 用 
于 中 国 地 磁 台 站 的 效果 进行 评估 , 提出 了 一 种 改进 的 
地 磁 Kp 指数 现 报 模 式 , 新 模式 能 有 效 识 别 地 磁 静 日 
变化 的 逐日 变化 特性 , 反映 地 磁 扰动 的 季节 效应 和 地 
方 时 效应 , 从 而 提升 了 Kp 指数 现 报 的 准确 性 . 通过 
XT EG, 检验 了 新 模式 的 效果 , 进而 分 析 了 Kp 指数 现 
报 的 误差 及 可 能 的 原因 . 


1 Kp 指数 现 报 模式 及 应 用 效果 


1.1 Kp 指数 现 报 模式 

Kp 指数 是 由 全 球 13 个 台 站 的 KK 指数 推算 得 
到 的 , 计算 过 程 包括 5 个 步骤 : (1) 提取 单 台 站 地 磁 
三 分 量 的 规则 日 变化 ; (2) 从 地 磁 实 测 值 中 去 除 规则 
日 变化 , 计算 单 台 站 的 3n 扰 幅 ; (3) 根据 扰 幅 与 KK 
指数 对 应 关系 , 计算 单 台 站 K 指数 , 0-9 共 10 
级 ; (4) 根据 K-Ks 转换 表 计 算 单 台 站 Ks 指数 ; (5) 
所 有 Kp 台 站 (分 布 如 图 1 Br) 的 Ks 求 平均 得 
到 Kp 指数 . 由 于 不 同 磁 纬 度 的 台 站 、 同 一 台 站 不 同 
季节 (使 用 劳 埃 德 天 文 季节 划分 法 , 一 年 划分 为 三 个 
季节 UU) 或 地 方 时 的 地 磁 扰 幅 平均 强度 均 不 同 , 在 历 
史 频 次 分 布 相同 的 原则 下 , 每 个 标准 台 站 都 制定 出 各 
自 扰 幅 与 K 的 对 应 关系 、K-Ks 转换 表 , 以 使 得 同 
一 时 段 各 台 站 Ks 指数 具有 相同 的 权重 , 所 有 标准 台 
站 的 Ks 平均 值 近似 反映 全 球 磁 情 721. 

Kp 指数 的 现 报 , 一 般 是 按照 Kp 指数 的 推导 
过 程 , 在 某 些 环节 上 采用 近似 假设 进行 计算 得 到 的 . 
在 计算 流程 上 , Takahashi 模式 与 Kp 指数 现 报 模式 
的 区 别 主 要 是 步骤 (1) 和 步骤 (4), 即 地 磁 规 则 日 变 
化 QDC (Quiet Day Curve) 的 提取 和 单 站 K 指数 
向 Kp 指数 的 转换 . 

地 磁 规 则 日 变化 QDC 的 提取 是 计算 Kp 指数 
最 重要 的 一 步 ， Menvielle 等 !7) 评估 了 多 种 计算 机 
标定 QDC 的 算法 , 并 与 手工 标定 结果 进行 对 比 , 发 
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图 1 Kp 指数 计算 使 用 的 134-835. SWPC 现 报 Kp 指数 使 用 的 9 个 台 站 及 北京 十 三 陵 地 磁 台 的 
位 置 分 布 (IGRF 2000 地 磁 坐 标 系 ) 
Fig.1 Locations of 13 geomagnetic field stations used for the derivation of Kp index, 9 stations used by SPWC for 


Kp index nowcast and the domestic station BMT (in IGRF 2000 geomagnetic coordinates) 


Ji FMI 方法 效果 最 好 . 但 是 FM 方法 无 法 直接 应 
用 在 Kp 指数 现 报 , 因为 在 求解 某 一 时 段 的 QDC 时 ， 
需要 使 用 该 时 自前 后 各 一 段 时 间 内 的 地 磁 实 测 值 , 这 
在 地 磁 指 数 实时 现 报 中 是 不 可 能 满足 的 . 为 解决 该 问 
题 , Takahashi 等 116) 经 过 分 析 评 估 , 采用 了 一 种 时 序 
ZME QDC 计算 方法 : 利用 当日 之 前 10 天 的 历史 
数据 , 对 全 天 0, 1, 2, …, 23 dE 24h, 求 之 前 10 天 
内 各 小 时 的 中 值 , 再 对 得 到 的 全 天 24 个 点 进行 傅 里 
叶 变换 , 滤 去 5 次 以 上 谐 波 , 插值 回 分 钟 精度 即 得 到 
当日 QDC. 相 比 Chapman 等 5] 提出 的 采用 当月 5 
个 磁 静 日 的 平均 值 计算 QDC 的 常规 方法 , Takahashi 
方法 是 一 种 改进 版 本 , 因为 其 不 再 依赖 于 对 磁 静 日 的 
人 为 判定 . 但 其 共同 假定 是 在 进行 时 序 色 加 的 时 间 
范围 内 QDC 的 形态 和 变 幅 无 明显 变化 . 

各 站 K 指数 向 Ks 指数 转换 是 为 了 消除 地 磁 
扰动 的 季节 和 地 方 时 差异 . Takahashi 现 报 模式 中 没 
有 区 分 季节 , 只 分 地 方 时 按照 线性 拟 合 关系 , 将 各 站 
的 K 指数 转 为 Ks. 

1.2 ”应 用 效果 

Takahashi 等 在 构造 Kp 指数 现 报 模式 时 , 使 用 
的 地 磁 台 站 的 地 磁 纬 度 均 在 43° 以 上 , 而 中 国 所 处 地 
磁 纬 度 较 低 , 空间 天 气 预 报 中 常用 的 北京 十 三 阶地 磁 


台 站 地 磁 纬 度 为 29"N, 直接 将 Takahashi 现 报 模式 
应 用 于 该 站 效果 不 理想 . 

K 1 列 出 了 利用 2000—2006 年 北京 十 三 陵 地 磁 
监测 数据 对 Takahashi 地 磁 Kp 指数 现 报 模 式 的 评 
佑 结果 . 根据 现 报效 率 定义 ， 


Pg-1- 


这 里 P HMR, O 为 实测 值 , O 为 实测 值 的 平均 值 . 
整体 上 分 析 , 现 报 值 与 实测 值 的 相关 系数 为 0.81, 现 
报效 率 为 0.70. 从 现 报 效果 的 地 方 时 和 季节 分 布 看 ， 
地 方 时 差异 非常 大 , 正午 附近 的 预报 效果 最 差 , 三 个 
季节 的 现 报效 果 有 一 定 差异 . 

为 提升 Kp 指数 现 报效 果 , 进一步 应 用 于 中 国 地 
磁 台 站 , 通过 改进 QDC 提取 方法 , 兼顾 地 磁 扰 动 的 
季节 和 地 方 时 分 布 , 提出 一 种 改进 的 Kp 现 报 模式 . 


2 改进 的 Kp 指数 现 报 模式 


2.1 地磁 规则 日 变化 QDC 实时 提取 方法 
地 磁 3h 扰 幅 是 下 类 指数 的 核心 变量 , 提取 扰 幅 
的 关键 是 从 地 磁 数 据 中 识别 非 磁 扰 成 分 . Bartels 等 
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表 1 Takahashi 现 报 模式 应 用 于 北京 十 三 陵 台 站 的 效果 
Table 1 Performance of the Takahashi Kp nowcast algorithm when 
applied to the BMT station 
E D 季节 季节 J 季节 全 年 

LC PE LC PE LC PE LC PE 
时 段 1 0.87 0.83 0.86 0.82 0.86 0.80 0.86 0.82 
时 段 2 0.75 0.71 0.81 0.70 0.74 0.61 0.77 0.68 
时 段 3 0.67 0.52 0.67 0.52 0.70 0.56 0.68 0.53 
时 段 4 0.69 0.48 0.68 0.47 0.76 0.55 0.71 0.50 
I BE 5 0.78 0.64 0.81 0.64 0.82 0.66 0.80 0.65 
时 段 6 0.87 0.80 0.87 0.79 0.88 0.78 0.87 0.79 
时 段 7 0.89 0.83 0.89 0.82 0.89 0.81 0.89 0.82 
时 段 8 0.90 0.83 0.91 0.84 0.87 0.79 0.89 0.82 
平均 值 0.80 0.70 0.81 0.70 0.82 0.69 0.81 0.70 


$È D 季节 为 1, 2, 11, 12 H; E 季节 为 3, 4, 9, 10 H; J 季节 为 5, 6, 7, 8 H. 时 段 1 到 时 段 8 分 别 对 应 北京 地 


方 时 02:00 LT 一 04:00 LT, 05:00 LT 一 07:00 LT, ---, 20:00 LT 一 22:00 LT, 23:00 LT—01:00 LT. LC 为 现 报 值 与 


实测 值 的 相关 系数 , PE 为 现 报效 率 . 


对 非 磁 扰 成 分 的 定义 是 非 K 变化 (non-K-varia- 
tions), 其 包含 了 太阳 静 日 变化 Sq 及 太阴 静 日 变化 L 
和 磁 扰 过 后 恢复 期 环 电流 等 磁 层 电流 缓慢 变化 Dia 
等 成 分 UU. Mayaud 等 给 出 了 识别 非 K 变化 的 一 
系列 规则 C, 符合 此 规则 的 磁场 变化 成 分 称 为 SR 
变化 , 即 规则 日 变化 . 本 文 沿用 Takahashi 模式 [79 
中 使 用 的 QDC (Quiet Day Curve) 表示 磁场 规则 日 
变化 曲线 , 将 Takahashi 模式 计算 QDC 的 方法 称 之 
为 Takahashi 方法 . 

QDC 具有 显著 的 逐日 变化 , 主要 归 因 于 太阳 日 
变化 Sq 的 逐日 变化 PO, S, 主要 取决 于 周 日 潮汐 引 
起 的 电离 层 发 电机 电流 , 其 具有 复杂 的 形成 机 制 , 但 
主导 变量 可 归结 为 电离 层 电导 率 和 潮汐 风 场 917221, 
电离 层 电导 率 及 潮汐 风 场 等 因素 的 逐日 变化 导致 5 
的 逐日 变化 P- 而 磁 层 电流 的 缓慢 变化 也 会 
对 QDC 引入 逐日 变化 的 成 分 P979... 图 2(a) 给 出 
的 是 2013 年 连续 两 个 磁 静 日 十 三 陵 台 站 H 分 量 随 
地 方 时 的 变化 , 表明 每 日 正午 附近 的 磁场 变化 均 非 常 
平缓 , 但 日 变化 幅度 和 相位 都 有 差异 . 

不 同方 法 识别 QDC 逐日 变化 的 效果 不 同 . 目前 
已 有 的 识别 QDC 的 方法 非常 多 , 按 基本 的 数学 方法 
可 分 为 统计 回归 法 、 频率 滤波 法 、 正 交 基 分 解法 等 几 


类 , IAGA 认可 了 USGS, AS, FMI 和 LRNS 等 4 种 
具体 方法 中 ,Menvielle 等 评估 的 结果 显示 , FMI 方 
法 与 人 工 标定 效果 最 接近 UT. FM 方法 是 一 种 变 窗 
口 回归 方法 : 在 计算 某 一 时 间 段 的 QDC 时 , 需要 衡 
量 这 段 时 间 前 后 的 磁场 变化 , 即 如 果 磁 扰 较 小 时 , 取 
较 小 的 时 段 内 的 磁场 平均 值 计 算 QDC; 而 当 磁 扰 较 
大 时 , 按照 磁 扰 规模 相应 扩大 时 段 计 算 QDCES-30. 
时 间 和 窗口 的 扩展 范围 按照 m 十 n 计算 . 其 中 : m 为 取 
决 于 地 方 时 的 常数 ; n — K33, 这 里 下 为 QDC 计算 
过 程 中 初步 得 到 的 K 指数 值 . 如 图 2(b) 所 示 , 由 于 
考虑 了 每 日 磁场 实际 变化 的 特点 , 该 方法 可 以 很 好 地 
区 分 QDC 的 逐日 变化 . 图 2(c) 表明 , 相 比 FMI 方 
法 , Takahashi 方法 无 法 识别 QDC 的 逐日 变化 . 

两 种 方法 可 以 结合 在 一 起 实时 识别 QDC. 对 于 
某 个 磁场 分 量 , FM 方法 计算 得 到 的 QDC 经 过 各 
个 时 段 时 间 窗 口内 实测 数据 的 平均 值 点 UT. mr 
QDC 本 身 , 计算 其 各 个 时 段 时 间 和 窗口 的 平均 值 然 后 
经 过 插值 、 渡 波 得 到 的 曲线 与 自身 重合 , 即 把 QDC 
代入 FMI 方法 将 得 到 QDC RA. 因此 , 如 果 先 利用 
过 去 数 天 数据 的 平均 值 初步 求 得 规则 日 变化 (下 文 
称 为 准 QDC), 再 把 准 QDC 中 地 方 时 与 当前 时 刻 之 
后 时 段 相对 应 的 值 取 出 , 连接 在 实测 数据 之 后 , 代入 
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图 2 20134 10 H 18 日 和 19 日 十 三 陵 台 站 地 磁 H 分 量 的 实测 值 (a), FMI 方法 的 QDC HAE (b) 以 及 
Takahashi 方法 的 QDC 计算 值 (c) 随时 间 的 变化 
Fig.2 Time variation of BMT's magnetic H component's observations (a) and QDC curves calculated by FMI (b) 
and Takahashi (c) method from 18 to 19 in October 2013 


FMI 方 法 , 则 在 实测 时 间 段 , 将 得 到 与 FMI 相同 的 结 
果 , 在 未 来 时 间 段 , 将 得 到 与 准 QDC 相同 的 结果 , 在 
两 部 分 数据 的 连接 处 , 计算 QDC 的 效果 受到 磁 扰 程 
度 和 消除 两 端 数 据 跃 变 方法 的 影响 . 在 磁 静 时 , QDC 
由 Sa 变化 决定 . Xu 等 P9 认为 , Sa 的 逐日 变化 主要 
表现 在 正午 附近 的 变化 幅度 上 , 而 相位 的 逐日 变化 大 
多 在 半 小 时 之 内 . 因此 , E QDC 在 数据 间断 点 前 一 
小 时 和 后 一 小 时 的 均值 变化 应 与 实测 数据 相应 的 变 
化 相差 不 大 , 变化 量 值 的 差异 取决 于 QDC 的 幅度 . 
因此 把 准 QDC 的 未 来 部 分 整体 上 下 平移 使 得 合成 
数据 在 间断 点 前 后 时 均值 的 变化 其 与 准 QDC 的 相 
应 变化 量 相同 , 即 可 在 如 下 意义 上 消除 间断 , 即 最 终 
得 到 的 QDC 与 准 QDC 在 两 区 间 连 接 处 的 变化 率 
相同 . 在 磁 扰 时 , 间断 点 前 一 小 时 的 磁场 强度 可 能 整 
体 变化 , 应 对 时 间 窗 口 进行 扩展 . 由 于 不 同 磁 纬度 台 
站 QDC 变 幅 与 地 磁 扰 动 幅度 的 比例 不 同 , 最 佳 扩展 
长 度 因 台 站 纬度 而 蜡 . 


图 3 给 出 的 是 利用 新 方法 生成 QDC 的 过 程 , 假 
设 当 前 时 刻 为 15:00LIB (北京 时 间 ), 所 现 报 3h 区 
间 为 12:00LTB 一 15:00LTB. 任 一 地 磁 分 量 实测 值 截 
止 到 15:00 LTB( 3a); 根据 过 去 数 天 实测 值 计 算得 
到 准 QDC( 图 3b); 对 准 QDC 进行 整体 平移 , 使 其 
与 实测 值 在 14:00 LTB—15:00 LTB 时 段 的 平均 值 相 
等 , 截取 平移 后 QDC 在 15:00 LTB 时 刻 以 后 的 部 分 ， 
计算 实测 值 在 14:00 LTB—15:00 LTB 时 段 的 平均 值 
时 , 应 根据 初步 得 到 的 C 指数 值 判断 向 前 扩展 时 间 
窗口 的 长 度 (图 3c); 因 FMI 方法 使 用 3 天 数据 计算 
中 间 一 天 的 QDC, 所 以 取 1.5 天 的 实测 值 , 1.5 天 经 
过 平移 的 准 QDC, 代入 FMI 方法 , 计算 得 到 中 间 一 
天 的 QDC, 截取 其 中 间 时 段 12:00 LTB—15:00 LTB 
的 部 分 作为 现 报时 段 的 QDC (图 3d). 图 4 给 出 的 
是 利用 新 方法 计算 QDC 得 到 的 结果 . 因 计 算 过 程 
按 3h 区 间 推 移 , 最 终结 果 由 所 有 3h 区 间 的 QDC 
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图 3 利用 新 方法 生成 QDC 的 过 程 . 15:00 BLT 为 当前 时 刻 , 12:00 LTB—15:00 LTB 为 所 现 报 的 3h KH. (a) 实 线 表示 
某 一 地 磁 分 量 实测 值 , 加 号 表示 每 小 时 的 平均 值 . (b) 虚线 表示 利用 过 去 数 天 观测 值 计算 的 准 QDC, 方 框 表示 每 小 时 的 
平均 值 . (c) 整体 移动 准 QDC 使 得 14:00 LTB—15:00 LTB 时 段 其 平均 值 与 实测 值 平 均值 相等 , 保留 移动 后 准 QDC 
15:00 LTB 时 刻 以 后 的 部 分 . (d) 将 实测 数据 与 移动 后 准 QDC 代入 FMI 方法 计算 得 到 24h 的 QDC, 
截取 其 中 12:00 LTB—15:00 LTB 时 段 
Fig.3 Example of calculating QDC by the new method, supposing that the present local time is 15:00 LTB and the 


three-hour-interval for nowcast is 12:00 LTB—15:00 LTB. (a) Solid line indicates observation of a geomagnetic field 
component, plus sign indicates hourly average of the observation. (b) Dashed line indicates quasi-QDC calculated 
from former several days’ observations, square sign indicates the hourly average of the quasi-QDC. (c) Moving 
quasi-QDC as a whole to a position where hourly average of the moved quasi-QDC is equal to that of the observation 
in 14:00 LTB—15:00 LTB. (d) Applying the FMI method to 3-day data, including 1.5 days of the observations before 
15:00 LTB and 1.5 days of the moved quasi-QDC after 15:00 LTB and getting the QDC of the central day. Taking the 
results corresponding to the interval of 12:00 LTB—15:00 LTB as the final QDC needed for the nowcast 


aë 2.2 K 指数 转换 为 Ks 指数 的 多 项 式 拟 合 
2436 [ — eio One 采用 多 年 样本 , 按照 上 述 方案 得 到 3h elit, 并 
sael MOG | 按照 各 台 站 转换 关系 得 到 KC 指数 估算 值 Ke JA, BU 
| 可 用 最 小 二 乘 回 归 方法 , 得 到 从 Kos 到 单 台 Kp fi 
S | 算 值 Kpa 的 转换 函数 , 用 于 未 来 现 报 . 在 不 同 地 方 
2353 时 及 不 同 季节 下 需要 统计 各 自 的 转换 函数 , 以 完整 将 
z 代 K-Ks 转换 表 . Takahashi 模式 采用 线性 拟 合 方式 ， 
a | | 区 分 全 天 8 个 3h 段 , 但 并 未 区 分 全 年 3 个 劳 埃 德 季 
18 r 20 节 . 新 模式 按照 式 (1) 使 用 多 项 式 拟 合 , 分 别 得 到 3 
a n 个 劳 埃 德 季节 下 8 个 3h 区 间 的 24 组 转换 函数 ( 见 


t áz — e E mpi 
Fig.4 Derived QDC using the new method (solid line), 较 大 , 不 同 季节 的 转换 函数 存在 一 定 差异 . 


compared with that of the FMI method (dotted line) Kp= f(Kes), Kpes = f (Kes). (1) 


X RAS 一 种 改进 的 Kp 指数 现 报 模式 
3 新 模式 现 报效 果 检 验 


3.1 “模式 对 比 

使 用 与 Takahashi 现 报 模式 相同 的 测试 数据 
(Fredericksburg 台 站 (FRD) 1998 年 2—10 月 数据 ) 
进行 测试 , 对 新 模式 与 Takahashi 模式 进行 对 比 , 结 
果 如 图 6 所 示 (图 6a 取 自 文献 16])， 从 图 6 可 以 
看 到 , 地 方 时 正午 附近 的 相关 系数 明显 提高 , 从 原来 
的 0.78 提升 至 0.88, 现 报 效率 从 0.62 提升 至 0.77; 
新 旧 两 种 模式 在 地 方 时 子夜 附近 现 报 值 的 相关 系数 
均 高 于 0.9, 新 模式 效果 略 有 提高 , 相关 系数 从 0.93 
提升 至 0.94, 现 报效 率 从 0.88 提升 至 0.89; 从 所 有 
样本 分 析 , 新 模式 Kp 现 报 值 的 离散 度 明 显 减 小 , 相 
比 旧 模式 , 相关 系数 从 0.87 提升 至 0.92, 现 报 效率 
从 0.77 提升 至 0.84. 可 以 看 出 , 新 模式 可 大 大 减 小 正 
F QDC 的 识别 误差 , 使 得 Kp 现 报 效果 明显 提高 . 

使 用 Takahashi 现 报 模式 9 个 台 站 ( 见 图 1) 中 8 
个 台 站 1998 年 的 数据 进行 多 台 站 测试 , 由 于 个 别 台 
站 数据 不 完整 或 有 错误 , 最 终结 果 取 各 台 站 结果 的 中 
值 . 相 比 Takahashi 模式 , 现 报 值 与 Kp 值 的 相关 系 
数 由 0.94 提升 至 0.96, 现 报效 率 由 0.88 提升 至 0.92. 
图 7 给 出 的 是 两 种 模式 的 最 终 现 报 值 与 Kp 值 的 相 
关系 数 随 世界 时 的 变化 , 可 以 看 到 , 新 模式 在 各 个 时 
段 均 优 于 Takahashi 模式 . 但 同时 两 种 模式 的 效果 都 
在 18:00UT 附近 较 差 , 这 可 能 是 由 于 所 使 用 的 台 站 
大 多 分 布 在 北美 ( 见 图 1), 正午 附近 的 现 报 偏差 没有 
完全 消除 . 图 8 给 出 了 两 种 模式 的 最 终 现 报 值 与 Kp 
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© 一 N w A ù ON N oo o 


159 


值 差异 的 频率 分 布 , 可 以 看 到 新 模式 有 8176 的 结果 
与 Kp 值 无 偏差 , 相 比 Takahashi 模式 的 7296, 提高 
了 996, 进一步 说 明 新 模式 具有 更 好 的 准确 性 . 两 种 
模式 均 有 约 99% 的 结果 与 Kp 值 偏差 不 大 于 1, 表 
明 多 台 联 合 现 报 可 相当 准确 地 反映 Kp 指数 量 值 . 
3.2 ”中 国 台 站 数据 的 现 报 效果 

利用 2000—2006 年 共 7 年 的 十 三 陵 台 站 地 磁 数 
据 进行 测试 , 分 别 评估 了 新 模式 与 Takahashi 模式 的 
效果 , 相 比 Takahashi 模式 , 新 模式 在 所 有 时 段 的 相 
关系 数 平均 值 由 0.81 提升 至 0.89, 现 报效 率 由 0.70 
提升 至 0.80. 图 9 给 出 了 两 种 模式 现 报 效果 的 比 
较 , 可 以 看 出 , 新 模式 相 比 Takahashi 模式 , 在 日 侧 时 
Et 05:00 LT—16:00 LT, 现 报效 果 提 升 较 明显 . 任 一 指 
标 均 存在 明显 的 地 方 时 变化 , 在 正午 时 段 11:00 LT— 
13:00 LT 效果 较 差 , 在 夜 侧 时 段 效果 较 好 . 图 9(c) 显 
示 , 均 方 根 误差 存在 季节 变化 , 在 季节 (春秋 分 ) 的 
正午 时 段 达到 最 大 值 , 现 报 结果 离散 度 最 大 . 


4 Wie 


图 10 给 出 了 使 用 2000—2006 年 北京 十 三 陵 台 
站 数据 , 采用 两 种 QDC 计算 方法 得 到 Kpes 均 方 根 
误差 随 Kp 的 变化 情况 ， 两 种 方法 的 变化 趋势 均 表 
现 出 Kp < 2 时 均 方 根 误差 较 高 ,Kp 为 2~7 时 均 方 
根 误差 递增 , Kp > 7 时 均 方 根 误差 减 小 . 新 方法 的 
效果 在 Kp < 6 时 优 于 Takahashi 方法 , 但 在 Kp >7 
时 基本 相当 . 


—— D14~17 


二 一 下 23 一 02 
一 JJ23 一 02 


图 5 Ke 到 Kpes 的 多 项 式 转换 函数 曲线 


Fig.5 Conversion function series from Kes to Kp, 
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图 6 使 用 FRD 台 站 1998 年 2—10 月 数据 时 Takahashi 模式 与 新 模式 的 现 报 结果 . (a) Takahashi 现 报 模式 给 出 
的 20:00 LT, 11:00 LT 和 全 天 现 报 值 与 Kp 值 的 对 比 ; (b) 新 模式 的 相应 结果 
Fig.6 Results are based on the data of FRD from February to October in 1998. (a) Distributions of real-time 
estimated Kp and Kp for 20:00 LT, 11:00 LT and all the local time given by Takahashi procedure. 


(b) Corresponding nowcast results given by the new procedure 


—— Takahashi procedure Cl Won OTO ETE 
---- new procedure 图 10 显示 出 两 种 QDC 计算 方法 用 于 中 国 台 站 
时 Kp 现 报效 果 的 差异 .新 方法 相 比 Takahashi Jy 
法 在 Kp = 6 时 有 明显 提升 , 但 在 更 强 磁 暴 时 提升 并 
不 明显 . Menvielle!” 指出 , FMI 方法 得 到 的 QDC 与 
实测 值 曲线 贴 合 太 紧 , 在 磁 扰 时 会 混入 规则 变化 以 外 
的 成 分 . Takahashi 方法 识别 的 QDC 虽然 与 实际 规 
则 日 变化 相 比 有 平均 10nT 的 偏差 59, 但 结果 较为 
E 稳定 . 新 方法 的 现 报效 果 取 决 于 FMI 方法 的 准确 性 ， 
- 因此 这 里 对 比分 析 了 FMI 方法 与 Takahashi 方法 的 
2 $ 8 9 5 2 效果 .图 11a) 给 出 了 两 种 方法 得 到 的 Kes 标准 关 
随 Kp 的 变化 , 在 Kp < 3 Bf, Takahashi 方法 的 结 
图 7 多 台 站 联合 现 报 结果 与 Kp 指数 相关 系数 果 较 差 , Kp > 3 后 Takahashi 方法 比 FMI 方法 的 结 
随 世 界 时 的 变化 果 稳 定 . 图 11(b) 给 出 了 两 种 方法 得 到 的 Kos 平均 
Fig.7 Variation of correlation coefficient between K pes 值 随 Kp 的 变化 , 显示 出 FMI 方法 得 到 的 结果 偏 小 ， 
and Kp using multi-stations network data 即 当 Kp 较 大 时 可 能 造成 漏 报 ， 但 对 于 最 终 现 报 结 
with Universal Time JR Kp 新 方法 相 比 Takahashi 方法 在 Kp « 6 时 有 
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图 8 多 台 联 合 现 报 结果 与 Kp 差 值 的 频率 统计 
Fig.8 Frequency distribution of the differences between Kp,, and Kp using multi-stations network data 
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图 9 2000—2006 年 十 三 陵 台 站 新 模式 与 Takahashi 模式 Kpe 及 Kp 在 不 同 季节 下 的 相关 系数 与 地 方 时 的 关系 (2), 现 报 
效率 与 地 方 时 的 关系 (b) 以 及 均 方 根 误差 与 地 方 时 的 关系 (c). 红色 为 新 模式 的 提升 部 分 
Fig.9 Local time variation of correlation coefficients between Kp,, and Kp (a), prediction efficiencies (b), and root 
mean square errors of the Takahashi procedure and the new procedure (c) using the BMT data from 2000 to 2006. 


Red bars indicate the improvement of the new procedure compared with the Takahashi procedure 
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更 好 的 效果 , 这 是 由 于 Kp < 3 时 Takahashi 方法 所 
得 结果 离散 度 较 高 而 造成 的 . 图 11(c) 给 出 的 是 Kp 
标准 差 随 两 种 方法 的 K。。 变 化 ， 相 比 FMI FA, 
Takahashi 方法 的 Kes 在 4~5 附近 时 Kp 的 离散 度 
RK. 图 11(d) 进一步 显示 ，Takahashi 方法 的 Kes 


1.4 


—— new method 


RMS(Kp.. — Kp) 


图 10 基于 2000—2006 年 十 三 陵 台 站 数据 所 得 新 方法 
和 Takahashi 方法 Kpe。 均 方 根 误差 与 Kp 的 关系 
Fig.10 Variation of the root mean square error of Kp. 
of the Takahashi method (dashed line) and the new 
method (solid line) with Kp using the data 
of BMT during 2000 to 2006 


--- Takahashi method 
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在 4-5 附近 时 Kp 的 平均 值 较 小 , 这 说 明 当 Kp f 
小 时 Takahashi 方法 更 可 能 出 现 虚报 , 使 得 对 中 小 
磁 暴 预警 的 准确 性 降低 ， 基 于 上 述 原因 ， 新 方法 相 
比 Takahashi 方法 在 Kp x 6 时 现 报 效果 更 好 , 但 
在 更 强 磁 暴 时 提升 不 明显 . 
4.2 Kp. 与 Kp 偏差 的 原因 

从 图 10 还 可 以 看 出 , 两 种 方法 的 现 报效 果 在 中 
强 磁 暴 (5 < Kp < T) 时 均 较 差 , 即 Kp。。 与 Kp 偏 
差 较 大 . 这 种 现 报 偏差 的 根本 来 源 是 3h 扰 幅 , 具体 
可 归结 为 两 方面 : 一 是 QDC 计算 方法 对 扰 幅 识别 的 
偏差 ; 二 是 相同 Kp 条 件 下 , 单 台 扰 幅 存 在 较 大 离散 
BE. 图 12 给 出 了 两 种 方法 现 报 结果 在 中 强 磁 暴 (5 < 
Kp < T) 下 标准 偏差 随地 方 时 的 变化 . 可 以 看 到 两 种 
方法 结果 的 离散 度 均 在 正午 时 段 11:00 LT—13:00 LT 
时 达到 最 大 , 在 17:00 LT—19:00 LT 时 有 一 个 局 部 极 
小 值 . 考虑 到 S4 变化 在 正午 附近 最 剧烈 ( 见 图 2), 从 
图 12 显示 的 现 报 结果 在 地 方 时 正午 附近 的 较 高 离散 
度 可 以 推测 , QDC 计算 方法 对 扰 幅 识别 的 偏差 是 中 
强 磁 暴 现 报效 果 较 差 的 重要 原因 . 此 外 , 地 磁 活动 存 
在 混沌 机 制 OY, 使 得 相同 Kp 条 件 下 , 单 台 扰 幅 本 身 
存在 离散 性 , 从 而 导致 现 报 结果 出 现 偏差 , 这 决定 了 
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图 11 利用 FRD 台 站 1998—2001 年 数据 得 到 的 Kp Kes 统计 结果 
Fig.11 Statistics on Kes and Kp by different QDC calculating methods using the data of FRD from 1998 to 2001 
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图 12 2000—2006 年 中 强 磁 暴 (5< Kp <7) 现 报 结果 
的 标准 差 随地 方 时 的 变化 
Fig.12 Local time variation of the standard deviation 
of the nowcast result by the new procedure (solid curve) 
and Takahashi procedure (dashed curve) during 
disturbed periods (5 < Kp < 7) based on 
the BMT data from 2000 to 2006 


单 台 现 报效 果 的 上 限 ， 如 图 12 所 示 , 两 种 方法 的 结 
果 在 17:00LT 一 07:00LT 同 有 0.75 左右 的 离散 度 . 确 
定性 混沌 所 致 扰 幅 的 离散 性 是 很 难 消除 的 ; 对 扰 幅 识 
别 偏差 的 消除 有 赖 于 进一步 改进 QDC 的 识别 方法 ， 
从 图 12 所 示 两 种 方法 现 报 结果 在 正午 附近 离散 度 的 
峰值 可 以 看 出 , QDC 的 识别 方法 仍 有 改进 余地 . 

另外 还 有 一 些 因素 可 能 会 对 现 报效 果 造 成 影响 ， 
例如 Kp 标准 台 站 的 时 区 分 布 不 均匀 , 主要 集中 在 欧 
洲 和 美洲 ( 见 图 1), 尽管 K-Ks 转换 表 则 在 消除 各 台 
站 磁 扰 强度 的 地 方 时 差异 , BARAR, Kp 指数 存 
在 明显 的 世界 时 变化 92: 磁 层 电流 系 的 空间 分 布 特 
征 , 单 台 站 的 平均 扰动 程度 存在 地 方 时 差异 . 

本 文 使 用 SuperMAG 指数 分 析 上 述 因 素 对 现 报 
效果 的 具体 影响 .SuperMAG 仿照 AE 指数 和 Dst 
指数 的 计算 方法 分 别 制定 了 极光 电 急流 指数 SME, 
部 分 环 电流 指数 SMR-00, SMR-06, SMR-12, SMR- 
18 和 环 电 流 指数 SMR (SMR-00, SMR-06, SMR-12, 
SMR-18 的 平均 值 ), 均 利 用 了 约 100 个 均匀 分 布 台 
站 的 数据 9633739, 反映 出 夜 侧 极 光 总 能 量 PI, AR 
同 地方 时 的 部 分 环 电流 强度 和 全 球 环 电流 的 平均 强 
BE 6 考虑 极 区 电流 系 ( 场 向 电流 及 相应 的 极 区 
电离 层 电流 ) 和 环 电流 这 两 个 影响 全 球 地 面 磁 场 扰 
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动 的 主要 电流 系统 , 使 用 SME, SMR-00, SMR-06, 
SMR-12, SMR-18 和 SMR 指数 分 别 作为 极 区 电流 强 
度 、 不 同 地方 时 的 部 分 环 电流 强度 和 对 称 环 电流 强 
度 的 衡量 指标 , 使 用 上 述 指数 的 3h 变化 幅度 分 别 表 
示 极 区 电流 及 不 同 地 方 时 的 部 分 环 电流 和 对 称 环 电 
流 的 3h 扰动 程度 . 
4.2.1 相同 Kp 条 件 下 极 区 电流 三 
环 电 流 强 度 随 世界 时 的 变化 

图 13 分 别 给 出 了 在 不 同 地 磁 条 件 下 与 Kp fü 
数 3h 区 间 对 应 的 时 间 段 极 区 电流 和 环 电流 的 扰动 
程度 随 世 界 时 的 变化 . 极 区 电流 的 扰动 程度 无 论 在 
磁 静 时 (Kp < 4) 还 是 中 强 磁 暴 时 (5 < Kp < 7) 都 
有 一 个 峰值 , 出 现在 09:00 UT—11:00 UT 时 段 内 . 环 
电流 的 扰动 程度 在 磁 静 时 随 世 界 时 变化 较 小 , 在 中 强 
磁 暴 时 存在 明显 的 两 个 峰值 , 分 别 位 于 12:00 UT 一 
14:00UT 和 18:00 UT 一 20:00 UT 时 段 . 考虑 到 Kp 
标准 台 站 有 约 一 半 集 中 分 布 在 欧洲 ，4 个 分 布 在 北 
美 ( 见 图 1), 从 图 13 的 结果 可 以 推断 Kp 标准 台 站 
在 正午 附近 容易 得 到 相 比 实际 全 球 磁 扰 强 度 偏 小 的 
结果 . 图 13(b) 的 两 个 峰值 一 定 程度 上 验证 了 这 一 观 
点 , 12:00 UT—14:00 UT 欧洲 台 站 处 于 正午 附近 , dE 
美 台 站 处 于 晨 点 附近 ; 18:00 UT 一 20:00 UT 时 欧洲 台 
站 处 于 和 昏 点 附近 , 北美 台 站 恰好 处 于 正午 附近 . 

因此 , Kp 台 站 的 经 度 不 均匀 分 布 导致 了 相同 Kp 
条 件 下 极 区 电流 与 环 电流 强度 随 世 界 时 的 变化 . DÀ 
为 极 区 电流 系 ( 场 向 电流 及 相应 的 极 区 电离 层 电流 ) 
和 环 电 流 是 两 个 影响 全 球 范围 的 地 面 磁场 扰动 的 主 
要 电流 系统 , 在 相同 Kp 量 值 下 全 球 范围 的 地 面 磁场 
扰动 强度 存在 世界 时 变化 . 
4.2.5 Kp 台 站 分 布 不 均匀 对 现 报 结果 的 影响 

42.1 节 的 结果 说 明 ，Kp 标准 台 站 的 时 区 分 布 
不 均匀 导致 了 相同 Kp 指数 下 全 球 范围 的 地 面 磁 拢 
强度 随 世 界 时 的 变化 .具体 到 单 台 站 ， 其 磁 扰 强度 
还 存在 随地 方 时 的 变化 . SMR-00, SMR-06, SMR- 
12, SMR-18 分 别 给 出 夜 侧 、 晨 侧 、 日 侧 、 展 侧 半 球 的 
环 电流 强度 , 显示 扰 时 不 同 地 方 时 的 环 电流 强度 BA. 
本 文采 用 十 三 陵 台 站 所 在 地 方 时 的 SMR 指数 的 3h 
扰 幅 来 大 致 衡量 相应 地 方 时 的 磁 扰 强度 . 08:00 LT— 
13:00LT 采用 SMR-12，14:00LT 一 19:00LT 3 
用 SMR-18, 20:00LT—01:00LT 采用 SMR-00, 
02:00 LT—07:00 LT 采用 SMR-06. 图 14(a) 为 不 
同 Kp 区 间 内 部 分 环 电流 指数 3 h 扰 幅 的 平均 值 随 
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图 13 2000—2006 年 中 强 磁 暴 时 (5 < Kp < T) 和 磁 静 时 (Kp < 4) SME 指数 (a) 和 SMR 指数 的 
3h WFE (b) 随 世 界 时 的 变化 
Fig.13 Universal time variation of the mean value of 3h range of SME index (b) and SMR index (a) for disturbed 
periods (5 < Kp < 7) and quiet periods (Kp < 4) from 2000 to 2006 
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Local Time (LT) 
图 14 2000—2006 年 SuperMAG 部 分 环 电流 指数 3h 扰 幅 平均 值 .(a) 暴 时 ( 实 线 ) 和 磁 静 时 (虚线 ) 随地 方 时 的 变化 ， 
(b) 暴 时 与 磁 静 时 的 比值 随地 方 时 的 变化 
Fig.14 Local time variation of the mean value of 3h range of SuperMAG partial ring current index from 
2000 to 2006. (a) For disturbed periods (solid curve) and quiet periods (dashed curve). 
(b) ratio of the disturbed periods and quiet periods 


十 三 陵 地 方 时 的 变化 ， 图 线 显 示 中 强 磁 暴 时 (5 « 识别 偏差 和 单 台 扰 幅 本 身 离散 性 的 影响 就 越 小 ;， 反 
Kp < T) 部 分 环 电 流 指数 扰 幅 存在 明显 的 地 方 时 变 — 之 ，QDC 识别 偏差 和 单 台 扰 幅 本 身 离 散 性 的 影响 
化 .对 各 个 地 方 时 ,中 强 磁 暴 时 (5 < Kp < 7) 和 就 被 放大 .， 图 14(b) 给 出 这 一 比值 随地 方 时 的 变 
HP (Kp < 4) 部 分 环 电流 指 数 扰 幅 的 比值 反映 出 ”化 , 其 变化 趋势 与 图 12 曲线 存在 反 相 关 性 , 特别 是 
此 地 方 时 附近 地 面 磁 扰 强 度 随 Kp 的 增长 速度 . E Æ 17:00LT 一 19:00LT 时 现 报 结果 离散 度 的 极 小 值 
值 越 大 扰 幅 分 布 的 范围 就 越 宽 ， 现 报效 果 受 QDC ”与 图 14(b) 的 峰值 相对 应 , 说 明 Kp 标准 台 站 的 不 均 
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匀 分 布 对 单 台 站 的 拢 时 现 报效 果 造 成 了 影响 , 
5 结语 


基于 提升 Kp 指数 现 报 效果 的 目的 , 提出 了 一 种 
改进 的 Kp 指数 现 报 模式 .新 模式 实时 识别 规则 日 
变化 的 方法 结合 了 Kp 计算 过 程 中 使 用 的 FMI 方 
法 与 SWPC 的 Kp 现 报 模式 中 使 用 的 Takahashi 方 
法 , 可 以 识别 规则 日 变化 的 逐日 变化 ; 新 模式 兼顾 了 
地 磁 活动 的 季节 效应 与 地 方 时 效应 , 使 用 24 组 多 项 
式 转换 函数 , 完整 替代 单 台 站 的 K-Ks 转换 表 . 从 整 
体 效 果 看 , 与 Takahashi 模式 相 比 , 新 模式 对 单 台 站 
现 报效 果 的 提升 非常 明显 ，Kp 台 站 FRD 的 相关 系 
数 从 0.87 提升 至 0.92, 现 报效 率 从 0.77 提升 至 0.84; 
中 国 BMT 台 站 的 相关 系数 从 0.81 提升 至 0.89, 现 报 
效率 从 0.70 提升 至 0.80. 新 模式 对 多 台 站 联合 现 报 
的 准确 性 也 有 所 提升 , 相关 系数 从 0.94 提升 至 0.96， 
现 报 效率 从 0.88 提升 至 0.92. 从 现 报效 果 的 地 方 时 
变化 看 , 在 地 方 时 正午 附近 新 模式 的 提升 效果 非常 显 
著 . 从 现 报效 果 随 磁 扰 程度 的 变化 看 , 新 模式 有 效 提 
升 了 中 小 磁 暴 (Kp < 6) 以 下 预报 的 准确 性 , 强 磁 暴 
的 预报 效果 与 Takahashi 模式 相当 . 

通过 评估 可 以 认为 , 条 件 有 限时 新 模式 可 直接 用 
于 中 国 单 台 站 实测 地 磁 数据 Kp 现 报 , 而 选取 时 区 跨 
度 较 大 的 多 台 站 数据 来 联合 现 报 能 够 获得 更 好 的 Kp 
现 报 效果 . 


致谢 所 使 用 的 Kp 数据 由 世界 数据 中 心 京都 地 磁 分 中 
心 提 供 , 使 用 的 SMR, SMR-00, SMR-06, SMR-12, SMR- 
18 指数 由 SuperMAG 提供 (http://supermag.jhuapl.edu). 
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